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We present results of numerical Self-Consistent Field (SCF) calculations for the equilibrium me- 
chanical unfolding of a globule formed by a single flexible polymer chain collapsed in a poor solvent. 
In accordance with earlier scaling theory and stochastic dynamics simulations findings we have 
identified three regimes of extensional deformation: (i) a linear response regime characterized by a 
weakly elongated (ellipsoidal) shape of the globule at small deformations; (ii) a tadpole structure 
with a globular "head" co-existing with a stretched "tail" at intermediate ranges of deformations 
and (iii) an uniformly stretched chain at strong extensions. The conformational transition from 
the tadpole to the stretched chain is accompanied by an abrupt unfolding of the depleted globu- 
lar head and a corresponding jump-wide drop in the intra-chain tension. The unfolding-refolding 
cycle demonstrates a hysteresis loop in the vicinity of the transition point. These three regimes of 
deformation, as well as the first-order like transition between the tadpole and the stretched chain 
conformations, can be experimentally observable provided that the number of monomer units in 
the chain is large and/or the solvent quality is sufficiently poor. For short chains, on the other 
hand, at moderately poor solvent strength conditions the unfolding transition is continuous. Upon 
an increase in the imposed end-to-end distance the extended globule retains a longitudinally uni- 
form shape at any degree of deformation. In all cases the system exhibits a negative extensional 
modulus in the intermediate range of deformations. We anticipate that predictions of patterns in 
force-deformation curves for polymer molecules in poor solvent can be observed in single molecule 
atomic spectroscopy experiments. 



I. INTRODUCTION 



The particular interest in the globular state of individual macromolecules (collapsed in poor 
solvent) and in the conformational collapse-to-swelling or unfolding transitions in individual poly- 
mer chains is motivated by the existing physical analogy between globules of synthetic polymers 
stabilized by solvophobic (attractive) interactions between the monomer units in poor solvents and 
the compact structures found for biopolymers, e.g., for globular proteins PQ. 

Despite this profound analogy and the fact that these objects share the same name, homopolymer 
globules differ markedly from protein ones. A globular protein, which is a copolymer with a large 
number of different amino acid residues, typically has a unique intra- molecular structure being a 
"aperiodic crystal" A globule of a flexible homopolymer can, in contrast, be better compared 
to a liquid droplet. In this case, the coil-globule transition is similar to the usual gas-liquid 
transition upon a decrease in the temperature (alternatively to the increase in the attraction 
between gas particles). The connectivity of the interacting monomers inside a chain ("linear 
memory") changes the characteristics of this transition to some extent; the transition keeps its 
phase nature but now becomes a continuous second order phase transition. The conformational 
and thermodynamic characteristics of a globule stabilized by monomer-monomer attractions are 
functions of the strength of this attraction, or in other words, are determined by the solvent quality. 
The density (segment concentration) of the globule grows monotonously upon a worsening of the 
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solvent quality from a very low value ~ TV" 1 / 2 in the Gaussian coil to a value of order unity in a 
densely packed globule. The remainder of the volume of a globule is occupied by solvent molecules. 
In the protein globule, on the other hand, the whole inner space is filled by the polypeptide chain 
and the globule is nearly dry. 

Recent developments in AFM force spectroscopy and optical tweezers techniques have made it 
possible to manipulate individual molecules and to subject them to mechanical force, for example, 
to stretch the chain by extending the distance between its two ends [3J HI [5J [7] . 

In essence, there are two possible scenarios for such extension. In the first case the governing 
parameter is the distance between two points in a macromoleculc (for example between the two end 
segments) which is fixed to a specified value or changed with a given velocity. The observable in 
this case is the reaction, or restoring, force. In the second case the applied force is fixed and plays 
the role of the control parameter. The observable in this case is the average end-to-end distance. 
In both scenarios a force-extension relation (or force-extension curve) is obtained which for finite 
chains not necessarily are identical. Note that in experiments with individual macromolecules 
the relevant distances are in the nanometer (nm) domain, whereas the force is in the piconewton 
(pN) range. In a recent review of Skvortsov et al. [8] it was demonstrated that the constant 
extension ensemble (first scenario) which is mostly used in experiments, leads to a remarkable 
reach deformation behaviour. 

There exists a large number of experimental works devoted to the extension of globular proteins. 
The objectives of such studies are to find, e.g., "weak spots" in a globule structure, or to discover 
a possible folding pathway. It has been shown in experiments on unfolding of proteins [S] that 
the force versus deformation curves may exhibit quite complex patterns and are essentially non- 
monotonic. 

The pioneering theory of unfolding homopolymer globules subjected to an extensional deforma- 
tion, was proposed by Halpcrin and Zhulina [10 . This theory envisions that a weak extensional 
deformation of a spherical globule, Figure [I] (a), produces a prolate ellipsoid of increasing asym- 
metry, Figure [I] (b) . The reaction (restoring) force was predicted to grow linearly at this stage. 
Under a moderate extensional deformation, a coexistence between a collapsed globular core and 
an extended "tail" takes place within a single macromolecule, Figure [l] (c) . The deformation in 
this "tadpole" regime is accompanied by a progressive unfolding of the globular core which occurs 
at an almost constant reaction force. This stage ends when a stretched string of thermal blobs is 
obtained, i.e. when the size of the tadpole's globular head becomes of the order of the thermal 
blob size, Figure [I] (d) . From this point onward the reaction force grows again with the following 
extension of the unfolded chain. 

Later on, Cooke and Williams [11] demonstrated the existence of a first-order conformational 
transition in the stretching of a collapsed (dry) polymer: at certain extension the chain suddenly 
unravels from the tadpole conformation, Figure[l](c), to the open chain conformation, Figure [T] (d) . 
In the force-extension curve this transition appears as a discontinuous drop in the force. The 
collapsed "head" in the transition point contains ~ TV 3 / 4 monomer units. The discontinuous drop 
in the force disappears, however, in the thermodynamic limit of infinitely long chain, N — > oo. 
Similar unraveling transition was found by Craig and Terentiev |12j who studied unfolding of 
globules made by semi-flexible polymers. 

The analogy between a homopolymer globule and a liquid droplet has been mentioned already. 
This analogy is also important for understanding the globule unfolding. The chain stretching, 
that is, the progressive increase in the end-to-end distance of the chain, is similar to increasing 
the available volume in the liquid-gas transition. The role of the gas phase is played in this case 
by the extended tail drawn out of the globule. Gas-liquid phase coexistence in a certain range 
of volumes (in our case - in the certain range of the given end-to-end distances) or at a constant 
pressure (in our case - at a constant reaction force) is well-known |13j . In the case of the polymer 
globule two phases coexist in a single macromolecule. This behavior is closely related to the 
Rayleigh instability in a liquid droplet [14] . In the case of a polymer in poor solvent an additional 
connectivity constraint of the solvophobic monomers in the chain |15j comes into play. Similar 
force-deformation patterns have been predicted for the unfolding transition in globular structures 
of amphiphilic associating copolymers [16J. The topological complexity, that is, the grafting of 
polymer chains onto a surface (a brush), introduces even more interesting features in the scenario 
of the unfolding transition [17], [15] , 

The aim of the present study is to develop a general theory of the polymer globule deformation as 
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FIG. 1: Unperturbed spherical (a) and ellipsoidal (b) globules, tadpole (c), and uniformly stretched (d) 
conformations. 

a function of the degree of polymerization N (number of monomer units in the chain) and the sol- 
vent quality. These two parameters determine the properties of the free unperturbed globule. Our 
work includes several stages. In the first part presented in this paper, we use the Scheutjens-Fleer 
self-consistent field (SF-SCF) lattice approach in its two-gradient version to study the evolution of 
the conformational and thermodynamic properties of a polymer globule upon a uniaxial extension 
(i.e. with an increase in the distance between the chain ends). We consider a wide range of N 
and polymer-solvent interaction parameter values. In the second part, an analytical Flory-type 
theory for the unfolding of a globule is advanced. The analytical theory uses, as input parameters, 
the properties of an unperturbed globule found with the aid of SF-SCF approach. This theory 
allows us to go beyond the limits of the SF-SCF modeling (first of all, to consider the limit of very 
long chains) and to calculate phase diagrams of the deformed globule as well as the properties of 
the globule in the transition point. These issues related to the phase diagrams, however, will be 
presented elsewhere. 

The remainder of the paper is organized as follows. In the section "Model and method" we 
introduce the SF-SCF approach. The results of the calculations are summarized in the section 
"Results" . In the section "Blob picture of globule deformation" we discuss the obtained results in 
terms of the classical scaling theory by Halperin and Zhulina for the extension of a flexible linear 
polymer chain collapsed in a poor solvent, and this is followed by the Conclusions. 



At the basis of the SCF approach is a mean-field free energy which is expressed as a functional 
of the volume fraction profiles (normalized concentrations) and the self-consistent field potentials. 
The minimization of this free energy leads for polymer chains to the Edwards diffusion differential 
equation |19j . which for an arbitrary coordinate system may be expressed as 



II. MODEL AND METHOD 



A. Theoretical background of the self-consistent field (SCF) approach 




G(r',r;n) 



(1) 
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The Green's function G((r', r; n) used in Eq. ([IJ is the statistical weight of a probe chain with the 
length n having its ends fixed in the points r' and r. The self-consistent potential u(r) represents 
the surrounding of the chain and serves as an external field used in the Boltzmann equation to find 
the statistical weight for each chain conformation. Consequently, the Green's functions G(r',r; n) 
that obey Eq. is related to the volume fraction profile of the polymer by a composition law: 

V{ ) ~ Er'Er" G(r",r';iV) " (2) 

In our case the ends of the chain are "pinned" at the points r' and r", so, the usual summation 
over the end points position should be omitted: 

The boundary conditions and incompressibility condition: <p(r) + tpsi?) = 1, where <ps(v) is the 
volume fraction of the monomeric solvent, provide constraints on the spatial solutions. The po- 
tential u(r) is local (i.e. there are no long-range forces) and depends on the local volume fraction 
tp(r) 

u(r) = u[<p(r)}, (4) 

The explicit concentration dependence will be specified below (see Eq. ([7])). Equations ([IJ, (J3j> , and 
Q make up the system of self-consistent field equations which is solved iteratively: one assumes 
an initial volume fraction profile ip(v), then computes the potential u(r) using Eq. Q, the set 
of the Green's functions, Eq. (JT|), and derives a new volume fraction profile ^'(r), Eq. The 
procedure is then repeated until the sequence of approximations ^(r) — > ip'(r) — > . . . converges to 
a stable solution <p*(r) = ip(r). 

To solve the self-consistent field equations rigorously, it is necessary to introduce a numerical 
algorithm. Such numerical scheme invariably involves space discretization (i.e., the use of a lattice). 
Here we follow the method of Scheutjens and Fleer (SF-SCF) [20] , who used the segment size a 
as the lattice cell size. The lattice sites are organized in layers, each of these layers is referred 
to with a single coordinate r. Within a layer, a mean-field approximation is applied, i.e., the 
volume fractions of the monomeric components and the self-consistent potential within the layer 
are constant. The way the sites are organized in layers depends on the symmetry in the system and 
must be preassumed. The approach allows for volume fraction and self-consistent field gradients 
between these layers. 

In order to consider the stretching of a single polymer chain (globule), it is necessary to use a 
two-gradient version of the SCF algorithm |21j , taking into account the symmetry of the problem. 
The natural geometry for this is a cylindrical coordinate system for which r = (r, z), Figure [2] 
In this case, all volume fraction profiles as well as other thermodynamic values depend explicitly 
on the radial coordinate r and the axial coordinate z. The system is rotationally invariant with 
respect to the z-axis and the mean-field approximation is applied along the angular coordinate. 
Therefore a lattice layer r = (r, z) represents a piece of a tube (a ring) of thickness and height both 
equal to the lattice unit length a. 

A polymer chain is represented on the lattice as a freely-jointed chain walk. The law of such 
walk on the cylindrical lattice is specified by setting the a priori transition probabilities, which 
are the statistical weights of the steps from the layer with the coordinate (r, z) to a neighboring 
one (r + i ■ a, z + j ■ a), where i, j G {—1,0, 1}, denoted by Ay. These statistical weights account 
for possibilities of going to the nearest neighbor (if i = but j =^ 0, or i ^ but j = 0), next to 
nearest neighbor (if i ^ and j =^ 0) or staying in the same layer (if i = j = 0). Since the number 
of sites in the layer depends on r, Ay is r-dependent but has no z dependence. The transition 
probabilities Ajj(r) obey the internal balance equation L(r)\\^(r) = L(r + a)A_ij(r + a) where 
L(r) = 7r(2r — a) /a is the number of sites in the layer with the radial coordinate r. The probability 
to go from the layer with radial coordinate r to that with r + 1 should be proportional to the 
contact area between these layers Xij(r) ~ A(r), the latter is given by A{r) — 2irra. Similarly, 



FIG. 2: Cylindrical lattice used in two-gradient SCF calculations. 



A_i j(r) ~ A(r — a). Using this area we may write 

Ai ' o(r)=Ai a^)' A ^ o(r)=Ai ^z(^r' 

Al ' ±l(r) = A2 ^LW' A - 1 ' ±l(r) = A2 L(r) 

For the transition probabilities in z direction Ao l ±i(?') = Ai, the probability to stay in the same 
ring Ao.o(f) = Ao- The probabilities Ay(r) should obey the normalization condition 

£ mo = i- (6) 

i=— 1,0,1 j=— 1,0,1 

In particular, this gives for the "limiting" transition probabilities Ao + 4Ai + 4A2 = 1 (at zero 
curvature). The set of Ao, A 1; and A2 determines the character of the walk on the lattice, i.e., the 
chain entropy and rigidity. For instance, setting A2 = eliminates next-to nearest neighbor steps. 

In the lattice approach, the iterative procedure of solving the system of the SCF equations is 
implemented as follows. Once the initial guess for the volume fraction profile, (p(r), is set, the 
self-consistent segment potential u(r) is calculated as follows 

^=log(l- i/3 (r))-2xMr)>. (7) 

where \ is the Flory-Huggins parameter describing the polymer- solvent interaction. The angular 
brackets in ^ denote a local layer average over the nearest and next-to-nearest lattice layers 

(X(r)) = (X(r,z))= J2 E *ij(r)X(r + i-a,z + j-a). (8) 

i=- 1,0,1 j=— 1,0,1 

Here the set of transition probabilities {Ajj(r)} introduced above for specifying the lattice walk is 
also used as the set of weight coefficients at averagung. 

The Green's functions G(r',r;n) can be computed from the recurrence relation expressing the 
fact that a chain of n monomers can be obtained by adding a monomer to a chain of n— 1 monomers: 

G(r',r;n) = (G(r',r;n-l))G(r;l), (9) 

where G(r; 1) is the partition function of a monomer which is simply given by the Boltzman law: 
G(r ; 1) = exp[— u(r) /fcgT] . The angular brackets denote here the weighted sum over the neighbours 
of the layer r = (r, z) which is calculated alike the nearest and next to nearest neighbor average, 
Eq. ([§. 
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Note that the set of {Xi,j(r)} (or, equivalently, {Ao,Ai,A2}) may be in principle different for 
density average and the Green's function calculation. Indeed, one might want to exclude the 
"diagonal" or next-to-nearest neighbor steps and set A2 to zero for the walk (i.e. in the calculation of 
G(r', r; n)), but take into account next-to-nearest neighbour interactions and set the corresponding 
A2 nonzero. 

B. Implementation of the SCF approach 

Numerical calculations using the SF-SCF approach described above were implemented using 
sfbox software developed in Wageningen University |22j . It allows to perform efficient high- 
speed calculations, even on a personal computer, sfbox uses the same set of {Ao,Ai,A2} both 
for modelling the chain walk, Eq. and performing the neighbor average, Eq. (j8j. For our 
calculations A = Ai = A 2 = 1/9 were chosen. This means that both nearest neighbor and next- 
to-nearest neighbor layers are taking into account in performing site average, Eq. ^ as well as in 
modeling the chain trajectory (calculating Green's function), Eq. 

In the calculations, a cylindrical lattice with limited size is used, therefore a simulation box, i.e. 
the range for r and z: r € [a, r max ] and z € [a, z max ], as well as the boundary conditions should be 
properly set. The simulation box should not be too large to make the calculations low time- and 
memory consuming. On the other hand, it cannot be too small, in order to avoid the edge-effects. 
In all cases we used r max larger then the radius of the unperturbed globule Rq. The system size 
in the z-direction depends on the value of the extension. 

The macromolecule is placed symmetrically in the box so that its ends are pinned at the z-axis, 
rx = (a, zi), rjv = (a, zn), equidistant from the corresponding boundaries: z\ — a = z max — zjq > 
2Rq. The distance between the chain ends is D — z^ — Z\. 

We followed a specific protocol to study the unfolding of the globule. The first calculation (the 
first run of sfbox) is made for the minimum distance between the ends of the chain D/a = 1. Then 
both zn /a and, correspondingly, the size of the box z max /a direction are successively increased by 
unity and the following run is made, etc. The solution obtained at i-th step is used as initial guess 
for (i + l)-th step. This procedure is repeated up to strong stretching, D/a ~ N (in practice - up 
to D/a « AT/2). 

Then a second series of runs is made. This series corresponding to the refolding of the globule. 
It starts from D/a ~ N. Subsequently both zjy/a and z max /a are successively decreased by unity 
down to D/a = 1. Again, the solution obtained at the previous step is used as initial guess in the 
following run. For each run we obtain the free energy and detailed profiles of the volume fractions. 

As a result, two free energy dependences on the chain extension and two sets of teh force- 
deformation profiles, i.e. for the forward (unfolding) and backward (refolding) runs, respectively, 
are obtained for each given pair {N,x}- Note that in the described scheme of SF-SCF modeling 
only one chain end corresponding to the AT-th monomer is moved, while the first monomer remains 
pinned in the same point (z\ remains unchanged). 

There are the following parameters in our model. The number of monomers units is equal to 
N. Unoccupied lattice sites are taken by a monomeric solvent (incompressible system) and the 
corresponding Flory-Huggins parameters for the polymer-solvent interactions is \. 

III. RESULTS 

We have performed SCF calculations for several chain length N ranging from 200 to 1000 and 
a series of x values ranging from \ — 0.8 to 2.0. In the SCF calculations, the results of two types 
are obtained: (1) macroscopic (large-scale) thermodynamic properties, first and foremost the free 
energy, and (2) polymer and solvent volume fraction distribution profile (i.e. local properties). 

First we consider how the free energy of the globule behaves with a consecutive increase or 
decrease in the distance D between the ends of the chain (i.e. upon the mechanical unfolding or 
refolding, respectively, of the globule) and obtain the force-extension curves - force vs. deformation 
dependences - for different A~ and x values. Then, by considering the evolution of the density pro- 
files, we analyze what conformational changes occur in the globule upon deformation and correlate 
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FIG. 3: Free energy as function of chain extension for N = 200 and various values of %■ Inset: free energy 
for x = 1-4 in the vicinity of the transition point. 



these with the features of the force-extension curve. 



A. Free energy curves 

Figure [3] shows an example of the free energy F as a function of the extension D calculated for 
N = 200 and a series of x-values. Both unfolding and refolding branches are plotted in this figure 
and, as one can see, two branches coincide both at small and high stretching values D, whereas 
at moderate extensions they differ within some range of D-values provided the Flory parameter x 
is high enough. This means that at moderate deformations, there exist two local minima of the 
free energy. The position and the width of the region where these two minima coexist depends on 
the x-value. The free energy is a functional of the density distribution F = F[ip(r, z)\ and the two 
minima correspond to different states of the deformed globule, and there are two different (locally) 
equilibrium volume fraction profiles <pi(r, z) and <~pi(r, z). The state that has a lower free energy is 
stable (global minimum), the other state is metastable. At the point where the two minima have 
the same depth, which occurs when the unfolding and the refolding branches intersect, D = D tr , 
the system suffers (on the mean-field level) a first-order-like phase transition. 

However, the two minima are separated by a free energy barrier. For finite chain length this 
barrier has a finite height and in reality the transition is smooth; the system can fluctuate between 
local minima. However in SCF calculations, where thermal fluctuations are suppressed, the system 
stays in the metastable state until it reaches the spinodal point, where the metastable minimum 
disappear, and then jumps to the other (global) minimum. Hence, the mechanical unfolding- 
refolding cycle exhibits in a SCF calculation a hysteresis loop. 

At lower values of x (x = 0.8 and 1.0) the dependence of the free energy on deformation is 
different: unfolding and refolding branches of the free energy superimpose completely. 
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D/a 



FIG. 4: Equlibrium reaction force vs. extension curves for the globule with N = 200 at various values of 
X- Inset: Reaction force calculated for forward (unfolding) and backward (refolding) runs for \ — 1-4- 

B. Force-extension curves 

Once the dependence of the free energy vs. extension is known, the reaction, or restoring, force 
can be found by differentiating the free energy F with respect to end-to-end distance D. The 
value of the reaction force is / = dF/dD and it acts against the extension. Figure |4| shows force- 
extension curves obtained from the free energy dependences on deformation, Figure 131 by numerical 
differentiation of the equilibrium free energy. The true transition point (in thermodynamic sense) 
D tr is obtained from the condition Fi(D tr ) — F2(D tr ) whereas the metastable states are excluded 
from these consideration. The kink in the free energy curves at the transition point gives rise to a 
jump in the reaction force. 

If the metastable states are taken into account, i.e. unfolding and refolding free energy branches 
are different, one obtains the hysteresis loop on the force-extension curve. As an example one of 
such curves is shown for x = 1-4 in the inset of Figure |1J 

Inspection of the force-extension curves clearly reveals three different deformation regimes. (1) At 
small deformations the reaction force strongly increases with extension. Note that at very small 
extension, D/a ~ 1, the force is negative: Fixing the chain ends to very small distances imposes 
high local concentration of the monomer units that is larger than the average equilibrium concen- 
tration in the globule; this is penalized by an increase in the free energy and this causes a repulsion 
between the end-monomers. (2) After reaching a peak value the reaction force slightly decreases 
(quasi-plateau) in a wide range of D/a and then drops down at the transition point D = D tr , where 
the equilibrium free energy dependence has a kink, Figure [3] Neither the extensive plateau nor the 
jump in the force are observed in the case of \ = 0.8. This is consistent with observation that the 
unfolding and refolding branches are equivalent and that there was no kink / transition point in 
this case. Interestingly, at \ — 1-0 t ne unfolding and refolding free energy branches coincide too, 
but the corresponding force-extension curve shows a pronounced jump. This can be attributed to 
the narrow width of the region where the two free energy minima coexist; this region cannot be 
resolved within our lattice approach since the latter has the resolution AD /a > 1. (3) After the 
jump, at strong deformations, the reaction force starts to grow again. The force-extension depen- 
dence in this regime is universal and is independent of the solvent quality \- This free energy and 
hence the reaction force have a purely entropic origin. As is discussed in Appendix, for a lattice 
chain, an analytical expression for the force dependence on deformation exists for this branch (in 
the whole range of chain extensions). This analytical result is shown in Figure [1] by the dashed 
line. 
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At this stage we mention that the transition between regimes (1) and (2) occurs continuously. 
In summary, an increase in x leads to (i) an increase of the reaction force at fixed end-to-end 
distance; (ii) the broadening of the quasi-platcau (regime (2) ), and (iii) an increase in the value 
of the jump of the force at the transition point. 

C. Density profiles 

To correlate the observed deformation regimes with (possible) conformational changes in the 
globule, profiles of polymer volume fraction distribution in the deformed globule should be ana- 
lyzed. As it was mentioned above, the SF-SCF approach gives access to these polymer density 
distribution profiles ip(r,z). The volume fraction is a function of two variables that can best be 
presented in two-dimensional contour plot, Figure [5] 

Another interesting and illustrative quantity that characterizes the conformations of the de- 
formed globule is the "integral" profile - the axial distribution of monomer units (the number of 
monomer units per z axis unit length) 

n{z) = Y,L(r)v{r,z). (10) 

A set of n(z) profiles for different extensions can be conveniently represented as a 3-dimensional 
plot, Figure [6] 

Figures [5] and [6] show the evolution of tp(r, z) and n(z) upon the imposed deformation for \ = 1.4 
and TV = 200. This choice of x an d N represents the typical case where three regimes of globule 
deformation, Figure [4] are well distinguished. One can see that at small extensions (regime (1)) 
the globule changes its shape from spherically symmetric (at D/a = 1), to an asymmetric one 
similar to that of a prolate ellipsoid (at D/a = 10). When the extension grows, the prolate globule 
conformation becomes unstable and the globule splits into a dense "head" and a stretched "tail" 
coexisting in one macromolccule (at D/a = 20), thus acquiring a tadpole conformation. This 
corresponds to regime (2) on the force-extension curve. When the distance between the ends 
of the chain increases, a redistribution of monomers between two phases occurs: the tail length 
grows and the head size decreases, this corresponds to the range of extensions from approximately 
D/a = 20 to D/a = 70. In spite of the decrease in the size of the globular head, the density 
of the globular core (i.e., except of the density in the diffuse interfacial layer) remains virtually 
constant. Close inspection reveals that the number of monomer units per unit length in the tail 
weakly increases. This is in accordance with a very weak decay of the reaction force / in the 
plateau regime, Figure |4j because the number of monomers per unit length is the inverse of chain 
extension: / ~ dz/dn = l/n(z). 

At a certain extension (at D/a 71 = D tr /a, in the transition point) the globular head 
disappears. Now the chain gets completely unfolded and one enters regime (3) (Figure |5jat D/a 
71 . . The disappearance of the globular head in the transition point leads to (i) a gain in the surface 
energy (since the interface disappears), (ii) a gain in the conformational entropy (since the chain 
tension decreases) and (iii) a penalty in the volume interaction free energy (since the monomer units 
that constituted the globular head get exposed to the solvent). The mutual cancellation of these 
three contribution determines the transition point. At larger extensions the tadpole conformation 
is metastable and the globule completely unfolds. The size of the disappearing globule in the 
transition point is still quite large (at D/a — 71 the numbers of monomer units in the tail and in 
the head can be estimated from Figures [5] and [6] 80). 

The density profile in the strong stretching regime has the shape of a homogeneous sphero- 
cylinder with local maxima in the points where the chain end are fixed. 

D. Special case: \ = 0-8, N = 200 

The case x — 0-8, N = 200 calls for special attention. The force extension curve, Figure [4] 
has an untypical shape compared to those for larger x values: there is no quasi-plateau regime at 



10 



10 



10 



Ss. 



10 



10 



10 



D/a=l 








0.5 



FIG. 5: Polymer density profiles <p(r, z) for the globule with N 
extensions. 
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FIG. 6: Number of monomer units per z axis unit length, n(z), for the globule with N = 200 and % — 1-4 
at various chain extensions. 

intermediate extensions and the force jump is missing. In this case the unfolding and refolding 
branches of the free energy vs. deformation coincide, see Figure [3] 

The analysis of the density profiles and the axial monomer distributions, Figures [7] and [8] re- 
spectively, demonstrates that in the considered case we encounter an essentially different unfolding 
mechanism: upon stretching the globule becomes more asymmetric and the density of its core 
decreases. What is also important, this transformation occurs continuously. 

The reason of such a behavior can be explained as follows. The globule has an interfacial 
layer with a finite thickness which increases as \ decreases (as the solvent becomes better for the 
polymer). In the case of a short chain, the width of the interface becomes comparable to the size 
of the globule (radius). It appears that it is thermodynamically more preferable to reduce the 
density in the globular core rather than to form an extended tail keeping the (high) core density. 
Therefore, instead of a redistribution of the monomer units between a collapsed (globular) and a 
stretched phase, these are transferred from the core to the interfacial layer, as D increases. 



E. Effect of the polymerization degree on the globule deformation 



Let us consider in more detail the influence of the degree of polymerization N and polymer- 
solvent interaction parameter \ on the force-extension curves. Figure [9]presents the force-extension 
curves calculated for different values of N and x- I n order to compare the data for different A^ 
the forces were plotted vs the reduced extension D/(Na). Figure [9] shows that for all considered 
cases, except of x = 0-8, N = 200, three deformation regimes are observed: (i) a linear force 
growth with an increase in D / (Na) at small extensions and (ii) a quasi- plateau in a wide range 
of D/(Na) values ending by a sharp decrease in / followed by (iii) a universal force-extension 
dependence of the ideal chain (see Appendix |A| Eq. (A4)). The position of boundaries between 
the regimes and the dependences of the restoring force on the deformation in regimes (i) and (ii) 
are determined by A" and \i the force-extension relation in regime (iii) depends only on N. The 
X-dependence was already discussed for A" = 200, Figure [4] Analysis of Figure [9] shows that an 
increase in either N or x has a similar qualitative effect on the shape of the force-extension curve. 
With an increase in x or iV the quasi-plateau region broadens mostly due to a noticeable shift 
of its right-hand boundary to larger D/(Na) values and the weak displacement of its left-hand 
boundary to smaller D/(Na). Therefore, it can be concluded that an increase in the chain length 
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FIG. 7: Polymer density profiles (p(r, z) for the globule with N = 200 and \ = 0-8 at various chain 
extensions. 



favors micro-segregation within the deformed globule. In particular we see that, in contrast to the 
case of N = 200 considered above, at x = 0-8 the globule formed by a longer chain unfolds via 
the formation of a micro-segregated (tadpole) structure with a subsequent drop in the reaction 
force. The height and the decaying slope of quasi-plateau also changes: with an increase in N 
and/or \ it shifts up and flattens. Upon a progressive increase in N the plateau level approaches 
some asymptotic height controlled by the value of \, whereas for given value of N it monotonously 
increases as a function of \- 

Taking into account the analysis of the conformations for TV = 200 as discussed above, we can 
presume that the picture of the unfolding of the globule as predicted by the SF-SCF method retains 
the main general features such as (i) the three stages of unfolding (three deformation regimes), 
i.e., (cf Figure [lj the growth of the force at small and strong deformations and its weakly decaying 
"quasi-plateau" behavior in the intermediate micro-segregated regime; (ii) a continuous microphase 
segregation (ellipsoid-tadpole) transition and the drop of the force in the second transition point 
where the micro-segregated structure becomes unstable and the globule completely unfolds. When 
the microphase segregation takes place, the head of the tadpole has a clear prolate ellipsoidal and 
not a spherical shape. On the other hand our SCF results clearly indicate that at moderately poor 
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FIG. 9: Force extension curves at various values of N and X- 



solvent quality the small globule unfolds without intra-molecular microphase segregated state. 



IV. BLOB PICTURE OF GLOBULE DEFORMATION 



In this section we discuss the above described results obtained using SF-SCF approach in terms 
of the blob picture of globule deformation. 
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A. A blob picture of polymer globule 



Let us consider again a single polymer chain with a degree of polymerization N immersed into 
a poor monomeric solvent. The chain is assumed to be intrinsically flexible, that is, the statistical 
segment length is of the order of a monomer unit length a, the latter coincides with the chain 
thickness. The binary attractive short-range (van der Waals) interactions between the monomer 
units are described in terms of Flory-Huggins interaction parameter x > 0.5. It is well known that 
in a poor solvent the polymer chain collapses into a spherical globule with a density ip which is 
a function of x- Under moderately poor solvent strength conditions, that is when (% — 0.5) -C 1, 
the monomer unit volume fraction in the globule scales as r ~ (\ — 0.5). The globule can be 
envisioned as an array of closely packed Gaussian thermal blobs, Figure [T] a. The blob size, 
£t ~ T~ l a is determined by the correlation length of the density fluctuations inside the globule. 
The size of an unperturbed globule scales as 



The chain fragments within a thermal blob retain Gaussian statistics, the number of monomer 
units per thermal blob scales as gt — (£t/a) 2 . The free energy of attractive monomer-monomer 
interaction per blob is of the order of thermal energy fc^T. Therefore the free energy of the 
globule scales (in the main term or in so called "volume approximation") as Fgio^^y/kgT ~ 
—N/gt — —Na 2 /£% ~ ~Nt 2 . A lower order correction in N arises due to the excess free energy 
of the interface between the collapsed globule and the (poor) solvent: the monomer units that are 
localized close to the interface exhibit in average more unfavorable contacts with the solvent than 
monomer units in the interior of the globule. In the scaling approximation an excess free energy 
of the order of ~ fc^T is attributed to each thermal blob localized at the globule interface. Hence, 
the excess free energy of the interface scales as Fgiotuie.s/ksT ~ (Ro/£,t) 2 — (Nt 2 ) 2 / 3 . This excess 
interfacial free energy stabilizes the spherical shape of a free (non-deformed) globule. 

Thus, thermodynamic properties of the globule are determined by the number of blobs (i.e. by 
the total number of blobs and by the number of surface blobs), while the globule size is controlled 
by both the number and the size of the blobs which are the functions of %. 



Following the scaling approach of Halperin and Zhulina [10 we now consider a deformed polymer 
globule as these occur when the end-to-end distance D is imposed. The extensional deformation 
D > 2Rq of the globule leads (under the constraint of the conservation of its volume) to an 
increase in the area of the interface, thus giving rise to an additional free energy penalty. The 
excess surface area of a weakly deformed globule, Figure [l]b, grows as AS ~ (D — 2R ) 2 as long 
as (D — 2R ) < R and here subscript "0" refers to the unperturbed globule. Hence, at weak 
deformations the restoring force that develops upon the extension of the globule, grows linearly 
with deformation, 



R ~ (N/g t yl^ t ~ (N/t) 1 / 3 



a 



(11) 



B. Globule deformation 




and can be presented as a string of Gaussian elastic blobs [23] 



D ~ (N/g){[ 



(14) 



of the size £ < £t- We remark that here we do not consider the non- linear elasticity effects that 
may occur due to the finite chain extensibility of the chains. The latter becomes important at 
limiting extensions, D ~ Na, see Appendix [A} 
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As has been noted by Halperin and Zhulina, in the intermediate range of extensions, 2Ro <C 
D <C (N/g t )£ t , the surface area and the excess interfacial free energy of a homogeneously elongated 
(cylindrical or prolate ellipsoidal) globule scales as F g i b u ie,s ~ D 1 / 2 . This leads to a decrease of 
the restoring force as a function of the extension. The non-monotonic behavior may be identified 
as a van der Waals loop, suggesting an intra-molccular co-existence of an extended chain fragment 
("tail") with a depleted spherical globular core, Figure [l]c. According to Halperin and Zhulina the 
unfolding of the globule in the intermediate range of extensions occurs at fairly constant force, 



f/k B T ~ 1/& 



(15) 



that corresponds to a tail consisting of a string of clastic blobs with a size equal to that of thermal 
blobs, £ t . If finite size corrections are considered, the coexistence plateau f/k B T ~ l/£ t ~ D° 
is replaced by a region where the force / is a weakly decreasing function of the extension D. 
Indeed, consider the globule in the tadpole conformation, Figure [l] c. Let us assume that the 
head of the tadpole has a spherical shape (that is, in the scaling analysis we can neglect the 
asphericity of the head) and is composed of n monomers. Then, the radius of the head scales 
as Rgiobuie — (n/gt) 1 ' 3 ^ — {na 2 ^) 1 ^ 3 - The free energy of the tadpole accounts both for the 
contributions of the globule and the tail: F ta dpoie — Fgiobuie + F ta u ■ The free energy of the globule 
contains volume and surface terms 



Fgiobuie 

k B T 



9t 



S, 



globule 

IT 



2/3 



(16) 



The tail consists of N — n monomers and can be represented as a stretched string of thermal 
blobs. The tail length is D — 2R g i buie- The number of blobs in the tail scales as (N — n)/g t ~ 
(N — n)(a/£ t ) 2 . Therefore the tail length can also be expressed as D — 2R g i b u ie — (N — n)a 2 /£ t . 
The free energy of the tail comprises an clastic contribution which is proportional to the number 
of thermal blobs in the tail 



Fum (D - Rgiobuie) 2 (N - n)a 2 



k B T (N-n)a 2 
As a result, the free energy of the tadpole is 



"tadpole 

k B T 

The restoring force is calculated as 



na 

Y 



i; 2 



2\ 2 /3 ( N -2n)a 2 







(17) 



(18) 



ftadpoie _ 1 dFtadpole/dn ^ 1 

k B T k B T dD/dn £ t 



2/3 



,1/3 



1 



(I) 



1/3' 



(19) 



and we notice that as D increases, then both n and thus / decrease. 



C. Ellipsoid-to-tadpole transition 

The location of the ellipsoid to tadpole transition can be estimated using simple arguments. A 
weak deformation of the globule, Figure[T|a, into the ellipsoid, Figure[l]b, produces a restoring force 



(12 1. The phase segregation inside the weakly deformed globule starts when this force becomes 



equal or comparable to the plateau force, Eq. (151 . This gives the threshold extension 



D ~ 2(6 + Ro) ^ 2& + (iVa 2 6) 1/3 ] =s 2(iVa 2 6) 1/3 (20) 



which size is of the same order as that of the unperturbed globule. 
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D. Tadpole-to-open chain transition 

The next step is to consider the transition from a microphase segregated tadpole, Figure [T] c, 
to the open conformation, Figure [T]d. In the transition point the free energies of the tadpole and 
open conformation should be equal: F ta dpoie = F chain- The free energy of the tadpole is given by 



Eq. (18 1. The free energy of the extended (open) chain is calculated similarly to F ta u, Eq. (17 1 

■ D 2 {N-nfa 2 



chain 



k B T ~ Na 2 ~ gN ' ^ 

where we have use the approximation that in the vicinity of the transition point the tail length is 
much larger than the globule size and, therefore, the tail size can be considered as being approxi- 
mately equal to the overall chain extension 

9t Si 

Equating the free energies of the tadpole and the open chain, we obtain that the number of 
monomers in the head of the tadpole (globule) in the transition point scales as 

n ~ N V 4 (^J 12 (23) 

or 

3/4 



n 
9t 



(0 



(24) 



that is, the number of blobs in the "minimal globule" is determined by the number of blobs in the 
whole chain. Scaling dependence of n vs. N in the transition point with the chain length N was 
obtained by Cooke and Williams [IT], who considered the limiting case of a "dry" globule. 

As a result, the size of the head at the transition point Rgiobuie — (na^t) 1 / 3 ~ A rl / 4 (^a) 1 ' /2 
depends on N and is much larger than the thermal blob size £ t . The jump in the reaction force is 

A/ ftadpoie f chain ^ D D ^ 1 at— 1/4 f £t\ (Oz,\ 

k^f~ k^f ~ (JV - n)a 2 ~ No 2 ~ a ' \H J ' 1 ' 

We see that the jump in the force at the transition point should increase with an increase in \ an d 
decrease with an increase in N thus vanishing in the thermodynamic limit N — oo . 



E. Arbitrary d case 

The above scaling analysis can be generalized in a way that takes into account an arbitrary 
dimensionality of the system d. In addition to the spherical globule case, which is the main subject 
of the present paper, corresponding to d = 3, we consider a polymer "bridging brush" immersed 
into a poor solvent and pulled by an external force [18] . which is effectively a one-dimensional 
system, d = 1. Similarly, as a two-dimensional case, d = 2, we can propose a cylindrical brush 
(bottle-brush) whose arms are simultaneously and equally pulled out in the radial direction (out 
of main chain or the grafting line) . 

If the globular head contains n monomers, the size of the globule scales as 
Rgiobuie ~ £t (na? / ^f) 1 d . Hence, the surface area per globule (per chain) is 

Sgiobuie ~ (Rgiobuie /Ct)^ 1 £f ~ (no 2 /^t)^ d This will generalize the surface contribution to 



the free energy of the tadpole, the second term in Eq. (18), the other contributions to F ta d P oie as 
well as to F c hain remain unchanged. The relation between the tail length and n is D — ^Rgiobuie — 
[N-n)a 2 /tl t . 
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Taking this into account, the generalized expressions for the reaction force in the tadpole regime 

-1/<T 



( 19 1 reads 



ftaa 



1 



i - 



i 



8 



i 

T 



i - 



d-l 



-l/d' 



(26) 



Since with an increase in D, n obviously decreases. Eq. (26 1 shows that f tadpole is a (weakly) 
decreasing function of the extension D for d = 2 and 3, whereas for d — 1 it is constant giving a 
true plateau on the force-extension curve. 

For the location of the transition of the tadpole to the open chain the expression ( 23 1 for the 
number of monomers in the tadpole's head is modified as 



N 3+1 



and for the force jump, Eq. (25 1, we have 

a/ _ i / e, 



k B T & \Na 2 



6 ' \n) 



(27) 



(28) 



Remark that in the case d — 1 corresponding to the planar polymer brush we recover the scaling 
n ~ TV 1 / 2 obtained earlier for the size of the microphase emerging in the brush capable to undergo 
a first order phase transition [2H H5J |5SJ |57j . 



F. Comparison of SCF results and scaling dependences 

The numerical results of the SCF modeling are in good qualitative agreement with the scaling 
laws obtained in the framework of the blob picture. In more detail, the result of the SCF theory 
have confirmed the existence of three regimes of extension of the globule as predicted by the blob 
model, namely (i) the deformation (extension) of the globule as a whole, (ii) the coexistence of 
a globular and unfolded (extended) phases, (iii) the disappearance of the globule and the further 
extension of the unfolded chain. The SCF results have further pointed to the possible departure 
from this scheme in the case of short chains and moderate values of the solvent quality \- 

A qualitative agreement of results of the SCF calculations and the blob model is also found 
for the dependences of the shape of the force-extension curve on the solvent quality and chain 
length. The parameters of the blob model are: (i) the blob size which decreases with \ an d (ii) 
the total number of blobs in the chain, the latter grows as a function of N and x- According to 
the blob model, the formation of an extended tail starts when the extension slightly exceeds the 
unperturbed globule size D ~ (JVa 2 ^) 1 / 3 . Correspondingly, the value of the ratio D/(Na) at this 
transition point should decrease upon an increase in x and N, exactly as is observed in the SCF 
modeling. According to the blob model, the height of the quasi-plateau on the force-extension 
curve should increase as a function of x (due to the decrease in £ t ) and N, that is, upon the 
increase in the total number of blobs. The latter determines the magnitude of the correction term 



in Eq. (19 1. A decrease in the number of blobs in the globular phase, n/gt, leads to a decreasing 
reaction force. Also, the jump-like decomposition of the globule at large D/(Na), as predicted by 
the scaling analysis is in agreement with the SCF analysis. 

At this stage it is necessary to realize that the blob model is not strictly applicable for the values 
of x an d N used in our SF-SCF simulations. The notion of a Gaussian thermal blob itself only 
has a meaning when the number of the monomer units in the blob, g t , is large enough. This is 
only the case at relatively small deviation from the 0-point, r = (9 — T)/T <C 1. The range of 
^-values considered in the present work obviously does not satisfy this demand. Furthermore, the 
total number of thermal blobs in the chain is ~ Nt 2 and this number should be large to use scaling 
arguments. Hence, the scaling parameter of the blob model is N 1 / 2 t, that requires at small values 
of t much larger chain lengths N as compared to those used in the numerical analysis. In the 
SCF modeling the value of N used in the calculations was restricted from above by computational 
reasons. 
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Admittedly, the blob model includes some simplifying assumptions as well. In particular, it is 
assumed that in the tadpole conformation, the globular head has a spherical shape. However, as 
one can see from the density profiles, Figure [5j the head has a prolate (ellipsoidal) shape rather 
than a spherical one. As a consequence, in the prolate globule, there appears a reaction force 
tending to restore its unperturbed spherical shape. In the tadpole this force is balanced by the 
elastic force from the tail. The analytical model that takes into account the non-spherical shape of 
the tadpole head and properly accounts for the force equilibrium will be considered in a follow up 
publication. We note, however, that the "spherical head assumption" does not affect the scaling 
dependences obtained in this section. 

V. CONCLUSIONS 

We have performed detailed SCF calculations of the equilibrium unfolding of a globule formed 
by a flexible homopolymer chain collapsed in a poor solvent and subjected to an extensional 
deformation. More specifically, we consider the conformational characteristics adopted by a chain 
with imposed end-to-end distances in a wide range of (poor) solvent qualities (expressed in terms 
of Flory-Huggins solubility parameter \) and polymerization degrees N . The fluctuating restoring 
force (i.e., the clastic force) is calculated as a function of the end-to-end distance. These results 
are collected in force-deformation curves. 

In accordance to predictions of the scaling theory by Halperin and Zhulina, we have observed a 
sequence of intra-molecular conformational transitions that occur upon an increase of the deforma- 
tion. We have found that there is a linear response regime (found at small deformations). This is 
followed by an intra-molecular micro-phase segregation regime (found at intermediate extension): 
here a uniformly stretched segment of the chain (a "tail") co-exists with a collapsed globular do- 
main (a "core" ) . A progressive increase in the end-to-end distance of the chain is accompanied by 
a systematic depiction of the globular core and a re-partitioning of the monomer units into the 
stretched tail. We have found that the unfolding of the globular core occurs at a weakly decreas- 
ing (fairly constant) reaction force, whereas the entropic elasticity is predicted to be recovered at 
strong extensions. 

The general shape of the force-extension curves obtained using the SF-SCF approach appears 
to be unconventional, showing a more or less extended region with an anomalous dependence with 
df/dD < 0, i.e. with a negative extensional modulus. In fact, f(D) curves exhibit the van der 
Waals loop, which is usually an indicator of the instability of a system. However, if the role of the 
governing parameter is played by the end-to-end distance, the system cannot avoid these states and 
passes through them step by step thus undergoing a sequence of intra-molecular conformational 
transitions. 

We have found that the intra-molecular co-existence occurs only for sufficiently long polymers 
and at large values of the \ parameter. In other cases, that is, for relatively short chains (number of 
segments of the order of 10 2 ) and mild solvent conditions x ~ 1, a uniformly stretched conformation 
is retained and no microphase segregation is predicted to occur in the whole range of extensional 
deformations. Nevertheless the force-extension curve exhibits a region with df/dD < 0. Hence, we 
predict that the phase diagram of the system in N,D or D coordinates contains an one-phase 
and two- phase regions and thus exhibits a critical point. 

Furthermore, we have analyzed two conformational transitions: (i) the first one is from the 
weakly elongated (ellipsoidal) globule to the tadpole, which consists of a ellipsoidal globular core 
coexisting with a stretched tail and (ii) the second one is from the tadpole conformation to that of a 
uniformly stretched chain - the unraveling transition discovered by Cooke and Williams [TT] . The 
first transition, that is the formation of a stretched intra-molecular micro-phase occurs continuously, 
whereas the second transition occurs as the first order phase transition and the reaction force drops 
abruptly down at certain elongation threshold. 

We anticipate that the patterns predicted by our theory should manifest in force-deformation 
spectra that can be obtained, e.g., by means of single-molecule AFM spectroscopy on end-grafted 
polymers collapsed upon a decrease in the solvent strength. The latter can be achieved by variation 
in temperature or in pH for thermo- or pH sensitive polymers, respectively. The most straight- 
forward way for a synthetic implementation of this system is to graft polymer chains to a solid 



19 



surface, for example by means of radical polymerization initiated at the surface. Typical degrees of 
polymerization obtained in controlled radical polymerization approach values as high as 10 2 — 10 3 . 

It is interesting to point to other polymeric systems that behave similarly to the stretched globule 
in a sense that these relieve the stress caused by an external field by "throwing out" a part of the 
chain as a new phase, thus forming a microphase segregated state. This is for example the case 
for a Gaussian chain compressed between two pistons. This chain undergoes an abrupt transition 
from a confined coil state to an inhomogeneous flower-like conformation partially escaped from 
the gap [28[ 29]. In the transition (microphase coexistence) region this system exhibits a negative 
compressibility, i.e. the reaction force decreases with an increase in deformation (corresponding 
to a decrease in a distance between the pistons, or chain squeezing) whereas at small and large 
deformation it grows. 

In the presented paper all results concern the deformation of a globule in the constant extension 
ensemble. In the conjugate constant force ensemble, one should expect that the tadpole conforma- 
tion will be unstable and the globule will unfold jumpwise, as a transition going from the ellipsoid 
globule directly to the open chain. This is confirmed by the results of molecular dynamics [30 
and Monte Carlo [31] simulations. However, in the case of lattice models or when the formation 
of helical conformation is possible a multistep transition can be observed [32, 33 . 
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APPENDIX A: CONFORMATIONAL FREE ENERGY OF (STRONGLY) STRETCHED 

CHAINS ON A CYLINDRICAL LATTICE 

The aim of this Appendix is to derive the conformational free energy of a chain in a cylindrical 
lattice as a function of the end-to-end distance D. It can be easily calculated if we consider an 
ideal chain extended by an external force (external field), i.e. we temporarily "switch" from the 
fixed stretching ensemble studied in this paper to the fixed force ensemble, that is the common 
Ansatz in statistical physics, and then switch back to the fixed deformation ensemble studied in 
this work. 

Suppose that the chain walking on the cylindrical lattice is subjected to a force / directed along 
the z-axis. For the sake of comparison with numerical results, we set the transition probabilities 
for the random walk according to SF-SCF model. Namely, the probability to make a step either 
in r or in z direction is given by Ai, the probability of a step in "rz" -direction, i.e. simultaneously 
changing both r and z coordinates by ±1 is A2 whereas the rest, Ao = 1— 4Ai — 4A2 is the probability 
to change the angular coordinate ip. Then the statistical weight of a monomer (of a link) is given 



Acknowlegdement 



by 



w = 



(A + 2A 1 )e° + (Ai + 2\ 2 )e fa ' kBT + (X 1 + 2\ 2 )e-i a ' kBT 




(Al) 



Then the partition function of the chain is 




(A2) 



The logarithm of the partition function gives us the Gibbs free energy, 
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Once the partition function (the Gibbs free energy) is known, the extension D corresponding to 
the force / can be found straightforwardly: 



D 



-fc R T- 



dGj_ 
df 



Na- 



2(A 1 + 2A 2 )sinh(^) 



l + 2(Ai+2A 2 ) 



cosh 



(fcir) 



1 



(A4) 



Eq. (A4| gives the dependence of the reduced extension (or the degree of stretching) D/(Na) on 
the applied force. This equation 

and we see that the force is a universal function of the reduced extension D /{Na). 

In order to return to the fixed extension ensemble, one can use the standard relation between 
the free energies in /- and D- ensembles (Gibbs and Helmholz free energies, respectively) 

F D =G f + D-f (A5) 

this gives the conformational free energy of the chain having its ends fixed at the distance D: 
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In Eq. ( A6 ) F is expressed as a function of / but together with Eq. ( A4 ) it parametrically defines 
F as a function of D. 

In the weak deformation limit we obtain 



Fchaii 
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k B T 4(Ai + 2A 2 ) Na 2 
So, we see that the elastic free energy has the Gaussian form 
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with k = l/4(Ai + 2A 2 ). The corresponding reaction force is 
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(A8) 



(A9) 



gives 



In the strong deformation limit smh.(fa/k B T) — cosh(fa/k B T) ~ \ exp(fa/k B T) and Eq. (A4) 

2(2fc — 1) • D/(Na)~ 



f chain 
knT 



log 



1 - D/(Na) 



(A10) 



The reaction force asymptotically tends to infi nity as D approaches Na. 

Comparing analytical force-extension curve ( A4 ) with SF-SCF numerical results (Figures [4] and 
[9] analytical curves are shown by dotted lines), we see a perfect correspondence in the uniformly 
stretched chain regime. 
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